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Alejandro Lopez-Bezanilla 1 ^ P. Ganesh 2 , and Peter B. Littlewood 1 3 
1 Argonne National Laboratory, 9700 S. Cass Avenue, Lemont, Illinois, 60439, United States 
2 Center for Nanophase Materials Science, Oak Ridge National Laboratory, 

One Bethel Valley Road, Tennessee, 37831, United States and 
3 James Franck Institute, University of Chicago, Chicago, Illinois 60637, United States 

Correlated band theory is employed to investigate the magnetic and electronic properties of dif¬ 
ferent arrangements of oxygen di- and tri-vacancy clusters in SrTiOs. Hole and electron doping 
of oxygen deficient SrTiOs yields various degrees of magnetization as a result of the interaction 
between localized magnetic moments at the defected sites. Different kinds of Ti atomic orbital 
hybridization are described as a function of the doping level and defect geometry. We find that 
magnetism in SrTi03_5 is sensitive to the arrangement of neighbouring vacancy sites, charge carrier 
density, and vacancy-vacancy interaction. Permanent magnetic moments in the absence of vacancy 
doping electrons are observed. Our description of the charged clusters of oxygen vacancies widens 
the previous descriptions of mono and multi-vacancies and points out the importance of the con- 
troled formation at the atomic level of defects for the realization of transition metal oxide based 
devices with a desirable magnetic performance. 


INTRODUCTION 

SrTiC>3 (STO) is a wide band gap insulator com¬ 
posed of nominally electrically neutral TiC>2 and SrO 
planes alternately layered in the form of a cubic per- 
ovskite. Alone or in contact with other materials STO 
has played a prominent role in the unrelenting growing 
field of oxide superlattices and transition metal oxide 
(TMO) electronics[J. Oxygen vacancies (Ov) are perva¬ 
sive point-defects in TMOs, and their presence impacts 
numerous of the material physical properties such as elec¬ 
trical conductivity or magnetism. Far from representing 
a drawback in the integration of TMOs in electronic com¬ 
pounds, the atomic control of Ov clusters as well as their 
density and distribution may lead to purposefully ma¬ 
nipulation of the magnetic and electronic properties of a 
TMO-based material [5]. 

Perovskites can accommodate a high density of Ov 
while keeping pseudocubic structure. Ordered Ov su¬ 
perstructures with a nominal stoichiometry of ABO2.5 
(brownmillerite) are formed upon removal of a substan¬ 
tial fraction of O atoms. Several experimental studies 
have revealed the presence of Ov clusters in alkaline-earth 
titanates, demonstrating that oxygen vacancies in STO 
can gather in highly n-type doped layers forming ordered 
domains with enhanced ionic conductivity^. Atomic 
precision synthesis of SrTi03_5 thin films demonstrated 
that vacancy-engineered structures can be artificially cre¬ 
ated. Neutron powder diffraction combined with syn¬ 
chrotron X-ray diffraction techniques have been used to 
determine the structure of complex-oxide Ti-based struc¬ 
tures identifying vacancy ordering arrangements in bulk 
samples[4]. Switching on/off behavior in Ti02-based 
resistive-random-access-memories has been found to be 
driven by the formation and disruption of conducting fil¬ 
aments of Ov[5]. 

Experimental techniques to increase or deplete charge 


carrier density Ov-rich areas allow the emergence of mag¬ 
netic phases in STO based compounds [6]. With a view to 
implementing TMO compounds in electronics, a detailed 
analysis of the properties of the magnetic and conduct¬ 
ing states present at the bulk of oxygen-deficient STO is 
of the uppermost interest. In spite of a considerable lit¬ 
erature on oxygen monovacancies, a qualitative system¬ 
atic description of oxygen multi-vacancies (mOv) in STO 
under different doping conditions has been lacking. Al¬ 
though previous studies have dealt with the formation, 
migration and isolation of Ov, the electronic and mag¬ 
netic properties of clustered Ov is still far from being 
fully understood. In particular, the understanding of the 
Ov geometry and concentration in bulk STO needs to be 
investigated to correlate the magnetic properties to the 
atomistic characteristics of the defected material. 


In this paper we analyse the large magnetic and elec¬ 
tronic tunability of oxygen deficient STO mediated by 
different doping rates of clustered oxygen vacancies. We 
use correlated band theory to predict magnetic behaviour 
in mOv and develop a comprehensive F-center analy¬ 
sis that accounts for the varied defect possibilities in 
SrTiOs-^. Unlike O monovacancies which may exhibit 
localized magnetic moments for particular doping con¬ 
centrations or vacancy densities [7], for any doping den¬ 
sity we find that isolated mOv always exhibit intrinsic 
magnetic moments. Both neutral and charged mOv con¬ 
figurations are considered to account for the effect of dop¬ 
ing in spin alignment. We find that electrons are confined 
to the vacated region in between the Ti atoms, produc¬ 
ing isolated F-centers whose magnetic moments are very 
sensitive to both the geometry and the material doping 
rate. 
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COMPUTATIONAL METHODOLOGY 

Density-functional-theory based calculations were con¬ 
ducted using the projector augmented-wave method [8] 
and the PBE-GGA exchange-correlation functional [9]. 
To improve the description of the electrons occupying the 
d-orbitals of the Ti atoms at the vacant site, a Hubbard-U 
correction (GGA+U) as implemented in the VASP code 
(T0iiT2| was included. The rotationally invariant method 
by Dudarev et al. [13] with an effective U e //=U-J=4.0 
eV was applied to capture the strong correlations. The 
electronic wavefunctions were described using a plane- 
wave basis set with an energy cutoff of 400 eV. Atomic 
positions were fully relaxed in large supercells until resid¬ 
ual forces were lower than 0.02 eV/A. Extra electrons 
were introduced or removed and compensated with an 
equally uniform background charge of opposite sign. 

RESULTS AND DISCUSSION 
Oxygen mono-vacancy 

We begin by recalling the origin of possible magnetism 
in the F-center created upon the removal of a single O 
atom from cubic STO. When STO is treated at low oxy¬ 
gen partial pressures and high temperatures, Ov are cre¬ 
ated and balanced by the generation of mobile charges 
according to the coordination of the oxygen atoms. In the 
perovskite structure, O 2- anions bond to Sr 2+ cations at 
the cube vertex and Ti 4+ cations at the cube center. In a 
pristine unit cell a Ti cation is 6-fold coordinated with the 
anions forming an octahedron. The highest solid state 
occupied bands contain 18 electrons primarily in oxygen 
p-orbitals hybridized with both cation orbitals. The low¬ 
est unnoccupied bands are composed of Ti 3d t 2 g states 
yielding an experimental band gap of 3.2 eV. This picture 
changes when an O atom is removed from the structure 
and the crystal symmetry is lowered from cubic Pm3m 
to the C 4 t, point group. In order to minimize the con¬ 
figurational energy the geometry in the vicinity of the 
vacant site undergoes a strong distortion including a lo¬ 
cal expansion of the volume of ~1%, a tilting of the Oq 
octahedra and an increasing of the Ti-Ti distance upon 
the formation of a new er-bond. Since the O vacates the 
crystal as a neutral atom the two electrons that accepted 
from the neighbouring Ti atoms are left in the defected 
structure, conferring to the Ov a double donor character. 

This traditional description of the localized magnetic 
moments in oxygen-deficient STO is based on density 
functional theory simulations that considered small su¬ 
percells and therefore favored the correlation of the defect 
states. In a previous study[7 we reported an extended 
model to explain the origin and conditions under which 
magnetic moments sitting at a single Ov are expected to 
be observed in STO. Therein we demonstrated by means 


of correlated band theory that magnetism in oxygen- 
deficient STO is not determined uniquely by the pres¬ 
ence of single Ov but by the density of free charge carri¬ 
ers. The local electron correlations between neighbouring 
oxygen deficient sites, the overlap of their wave functions 
and the structural rotation of the perovskite octahedra 
determine the appearence of a magnetic charge distribu¬ 
tion confined in the region between two Ti atoms, giving 
rise to magnetic moments which are further coupled to 
itinerant conduction band (CB) d-electrons. We found 
that isolated Ov are nearly hydrogenic double donors, 
while at higher carrier density, the Ov distorts to accom¬ 
modate a single bound electron, and therefore a local 
magnetic moment. 

Side-by-side oxygen di-vacancy 



FIG. 1 . From a) to e), electronic band diagrams displaying 
the evolution of the electronic states with increasing level of 
hole doping in a supercell containing two side-by-side oxygen 
vacancies as schematically shown in the inset, g) The for¬ 
mation of two magnetic states localized at the vacancy sites 
with opposite spin alignment is robust against doping, h) One 
electron at a vacant site is preferred over sharing the defect 
wave function on both site for three-fold hole doping. 

Ov are known to tend to cluster in specific patterns de¬ 
pending on their location either on the bulk or by the slab 
surface [13], Specific Ov arrangements such as line of Ov 
are strongly favored and form preferentially along certain 
crystallographic directions |15j. This results in a chain of 
Ti 2+ cations laterally bonded to only two O atoms. Also 
a side-by-side ordering where Ti-Ti bonds align in paral- 
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lei were found to be the most stable magnetic configura¬ 
tion in Ref HI], where both types of arrangements were 
found very close in energy. Upon considering several ge¬ 
ometries of the supercell containing the Ov clusters we 
determined that 3x4x7 and 3x3x11 supercells were suf¬ 
ficiently large (depending on the defect geometry) for an 
accurate description of the different spin alignments at 
the Ov sites while accounting adequately for the density 
of states, the octahedra rotation and the electron count¬ 
ing. Due to the reduced lateral extension and the long 
tail of the F-center wave function along the direction of 
the vacancy ordering, the longer dimension of the su¬ 
percell is selected in the Ti-Ti axis. Due the geometry 
distortion the vacancies may be correlated through the 
elastic constant of the network, but the large size of the 
supercell guarantees that the wave functions of periodic 
images do not overlap either laterally nor in the Ti-Ti 
bond direction. 

To investigate systematically the effects of Ov con¬ 
centration on the magnetic properties of STO, we start 
by analysing a di-oxygen vacancy (di-Ov) configuration 
where two side-by-side in-plane O atoms were removed 
yielding two parallel Ti-Ti bonds, as depicted schemat¬ 
ically in the inset of Figure [T] The removal of two O 
atoms releases four electrons that are accomodated in 
the defected bulk similarly as explained above for the 
monovacancy. One electron is trapped at each of the 
two vacant sites upon hybridization of the 3d 2 2 and 4p z 
atomic orbitals of each Ti atom, while the remaining elec¬ 
trons delocalize at the CB. This results in two magnetic 
states with opposite spin alignments and the same energy 
of ~ -0.5 eV, as exhibited in the band diagram of Fig¬ 
ure The Ti-Ti distance decreases a 1.3% to 3.91 A. 
A real-space representation of the F-centers localized in 
between the Ti atoms with different magnetizations are 
depicted in Figure [l)y Aligning both magnetic moments 
with same orientation entails an increase of the formation 
energy of 27 meV. 

This picture is barely modified if the system is posi¬ 
tively doped by removing one electron. The electron is 
substracted from the CB, as shown in Figure[l]), and the 
vacancies remain negatively charged. Additional removal 
of a charge carrier from the CB leads to a semiconducting 
system with the magnetic moments intact. This result 
is similar to the single-hole doped mono-Ov [7j, and sug¬ 
gests that up to n-hole doping (n= Ov concentration) the 
defect levels are pinned at ~ 0.5 eV below the bottom of 
the conduction band. Further removal of a third electron 
localizes the remaining electron in a vacant site. It is 
worth noting that the charge sits solely as an F-center in 
between a pair of Ti cations separated by 3.93 A, whereas 
the other two Ti atoms move apart from each other a 
distance of 4.1 A and slightly contribute to the magnetic 
moments with a small hybridization of the 3d,2 and 3d xz 
orbitals (see Figure |TJi) . If both vancant sites are forced 
to share the single electron wave function, the single state 


observed in Figure[l|l splits yielding two states with equal 
energy and different spin if the sites couple antiferromag- 
neticly, or with different energy and same spin alignment 
if coupled ferromagneticaly. In both cases the formation 
energy of the defected system increases 200 meV. Finally, 
removing 4 electrons the oxygen deficient structure be¬ 
comes a wide gap insulator as bulk STO. On the contrary, 
negative doping injects charge carriers in the bottom of 
the Ti d-bands and consequently enhances the metallic 
character by shifting up the Fermi level higher in the CB. 

Head-to-head oxygen di-vacancy 

An alternative arrangement of two Ov is a head-to- 
head configuration, which causes the reduction of two 
Ti 4+ cations to two Ti 3+ cations bridged by a Ti 2+ 
cation: -Ti 3+ -Ov-Ti 2+ -Ov-Ti 3+ - (see inset of Figure |2j. 
This arrangement is very close in energy to the side- 
by-side configuration discussed above and the computa¬ 
tional approach becomes decisive in determining the real 
ground state |16j. Despite all that, Ov in the proximity 
of domain wall and grain boundaries have been identi¬ 
fied as energetically preferred when segregated in linear 
clusters [Ti]. Starting from the ideal structure of STO the 
creation of this type di-Ov maintains the original frame¬ 
work of the perovskite structure although some octahedra 
tiltings around the defective sites are observed. 

According to the symmetry and the oxidation rate, 
different types of atomic hybridization can be observed. 
The band diagram of Figure[2]i shows that three electrons 
are trapped in three in-gap states whereas the fourth elec¬ 
tron is delocalized in the CB, rendering the system metal¬ 
lic. The projected density of states (PDoS) onto the Ti 
atoms is employed to assess the degree of electron local¬ 
ization on each atom, giving the the contribution of each 
cation orbital to the total density of states at a particular 
energy. Figure [3] shows that the magnetic state at ~-1.7 
eV is a hybrid of the 4s and 3d 2 2 orbitals of uniquely the 
Ti 2+ cation. The same two orbitals along with the p 2 and 
the 3d 2 2 orbitals of the Ti 3+ cations define the localized 
state at ~-1.5 eV, which has opposite spin to the latter. 
Finally, the shallow state at ~-0.2 eV is the result of mix¬ 
ing the p 2 orbitals of both type of cations and the 3d 2 2 
orbital of the Ti 3+ cations. A real-space representation 
of the net charge occupying the space along the Ti-Ov- 
Ti-Ov-Ti axis is given in Figure^. The total magnetic 
moment of the defect is 1 /r b ■ The distance between the 
Ti 2+ and a Ti 3+ cation is of 3.90A, and from a Ti 3+ to 
the Ti 4+ cation along the defect axis is of 4.00 A. 

The removal of an electron out of the 2676 contained in 
the supercell has extraordinary consequences on both the 
local chemistry and the magnetism of the di-Ov. Indeed, 
unlike the oxygen monovacancies and the di-Ov presented 
above, instead of a CB electron a localized state is re¬ 
moved when an electron is extracted from the supercell 
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Neutral (2 xOv) 1+ (2 xOv) 2+ (2 xOv) 3+ (2 xOv) J+ (2 xOv) 1 ” 



FIG. 2. From a) to e), electronic band diagrams showing the 
evolution of the electronic states with increasing level of hole 
doping of a supercell containing two head-to-head oxygen va¬ 
cancies. f) corresponds to a negatively doped configuration. 
The metallic or insulating behaviour, the total magnetic mo¬ 
ment, and the distribution of the magnetic states strongly 
depend on the doping level. The real-space representation of 
the net charge density of the undoped (neutral) configuration 
is displayed in g), for singly hole doped in h), for two-fold hole 
doped in i), and for three- and four-fold hole doped in j). The 
latter also represents a hybrid 4s3d,2 state for a di-oxygen 
vacancy in a supercell with four electrons removed. The ex¬ 
cess of spin-up with respect to spin-down charge density in a 
supercell is plotted in g), h) and j) for an isosurface value of 
1.7 x 10 _3 e“/A 3 , and in i) of 1.7 xl(T 4 e _ /A 3 . 


as observed in Figure [2 Jd. Interestingly, instead of reduc¬ 
ing the magnetic moment removing an electron enhances 
the total magnetic moment up to 2 fiB- The explanation 
of such a drastic change relies on the hybridization un¬ 
dergone by the Ti cations’ orbitals at the defect sites. As 
shown by the PDoS in Figure [3j the Ti 3+ cations has no 
contribution to the two in-gap states and, in turns, the 
central Ti cation orbitals accomodate two electrons. The 
deeper in-gap state at ~-1.5 eV is formed by a hybridiza¬ 
tion of the 4s and 3d z 2-orbitals, whereas the orbital 3d xz 
accomodates by itself alone an electron at ~-0.9 eV. No¬ 
tice that the molecular-like orbital formed upon the com¬ 
bination of the three cations vanishes and only one cation 
holds the magnetic moment of the entire system (Figure 
§r). The distance between the central cation holding the 
excess of charge and the neighbouring Ti 3+ cations along 
the defect line increases to 4.01 A, while the Ti 3+ -Ti 4+ 
distance is reduced to 3.93 A. 

Further removal of an electron renders the system a 
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FIG. 3. Spin-resolved projected density of states (PDoS) on 
the 4s, 4p z , 3d„ z and 3d z 2 orbitals of the Ti 3+ and Ti 2+ 
cations at the defect shown in Figure [I] consisting in two 
oxygen vacancies aligned in a SrTiC >3 slab. From upper to 
lower panels, the PDoS correspond to the neutral, single, two-, 
three-, and four-fold positively doped configurations respec¬ 
tively. 


non-magnetic insulator. The charge is withdrawn from 
the CB and a magnetic state with similar spatial ex¬ 
tension to the neutral configuration state (Figure [2j) is 
formed by the concurrent addition of the 4s, 4p z , and 
3 d z 2-orbitals of both the Ti 3+ and Ti 2+ cations. As a 
result, two magnetic states of opposite spin sign lay in 
the middle of the band gap. Observing the Figure W 
one might conclude that an additional extraction of an 
electron from the supercell removes all localized in-gap 
states. But a detailed inspection of the corresponding 
PDoS diagram shows that a flat state in resonance with 
the valence band (VB) states and formed upon hybridiza¬ 
tion of the 4s and 3d z 2 central Ti cation orbitals is re¬ 
sponsible of the supercell magnetic moment of 1 f^B- At 
high energy, close to the CB, the antibonding counter¬ 
part is observed. This spatially localized state is plotted 
in Figure [2], where the geometry is barely modified with 
respect to the previous singly doped structure. 

Interestingly, this is the same state that is observed 
when a total of four electrons are removed from the sys¬ 
tem. Contrary to what observed in the side-by-side di- 
Ov, the removal of all the Ov doping electrons does not 
necessarily entails the absence of a localized state and 
a magnetic moment, but the formation of a permanent 
magnetic moment sitting on the central Ti 2+ cation. Fi¬ 
nally, negative doping of the supercell with an additional 
electron shifts the Fermi level enhancing the metalic char¬ 
acter of the oxygen deficient system, as shown in Figure 
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Head-to-head oxygen tri-vacancy 
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FIG. 4. From a) to e), electronic band diagrams showing 
the evolution of the electronic states with increasing level of 
positive doping of a supercell containing three head-to-head 
oxygen vacancies. The metallic or insulating behaviour, the 
total magnetic moment, and the distribution of the magnetic 
states strongly depend on the doping level. Panels from g) to 
j) display a real-space representation of the net charge density 
of the SrTi 03_5 configurations. Real-space representations 
of the net charge density of the SrTiC> 3_,5 configurations are 
in display in f) for neutral, single, and two-fold hole doped 
defected supercell, and in g) and h) for three- and four-fold 
hole doped configuration respectively. The excess of spin- 
up with respect to spin-down charge density in a supercell is 
plotted in f) for an isosurface value of 2.6 xl0 -o e - /A 3 , and 
in g) and h) of 1.6 xlO~ 3 e“/A 3 . 

Next, we follow the same procedure to analyse the 
changes in the magnetic moment of three aligned Ov 
upon doping the supercell. Removing three O atoms as 
depicted in the inset of Figure [4j the supercell is left 
with two Ti 3+ and two Ti 2+ cations. This type of defect 
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FIG. 5. Spin-resolved projected density of states (PDoS) on 
the 4s, 4p z , 3d„ z and 3d z 2 orbitals of the Ti 3+ and Ti 2+ 
cations at the defect shown in Figure [4] consisting in three 
oxygen vacancies aligned in a SrTiC >3 slab. From upper to 
lower panels, the PDoS correspond to the neutral, single, two-, 
three-, and four-fold positively doped configurations respec¬ 
tively. 


exhibits a strong doping character due to the six elec¬ 
trons that are released. Figure [4ji shows that the neutral 
configuration corresponds to a metallic system with two 
electrons in the CB and two groups of flat bands each 
containing two in-gap states energetically degenerate and 
with opposite spins. Thereby the supercell exhibits a zero 
overall magnetization with several magnetic moments lo¬ 
calized in the Ti cations. Figure [4^ displays the Ti 2+ 
cations holding electronic charge with opposite spin to 
the Ti 3+ cations, and Figure[5]the substantial differences 
in the contribution of the cation orbitals to the local¬ 
ized states. The PDoS reveal that the four in-gap states 
are formed by the combination of the 4s, 4p z , and 3d z 2- 
orbitals of both types of cations. The weight of the Ti 3+ 
p~-orbitals in the hybridization is larger than those of the 
Ti 2+ cations, and conversely for the Ti 2+ 3d, 2 -orbitals. 
The degeneracy in energy of the magnetic moments is 
lifted when the Fermi energy shifts down when one elec¬ 
tron is removed from the conduction band (Figure §>)■ 
The defected system becomes insulating when a second 
electron is removed. Note that in no case the distribution 
of the net charge density is modified by the hole doping 
and both configurations exhibit the same magnetic mo¬ 
ments as the neutral configuration, and an overall zero 
magnetic moment. Also the distances between cations 
along the defected line are conserved through the three 
doping strengths: The distance between the Ti 2+ cations 
increases to 4.30A, from a Ti 2+ to a Ti 3+ cation reduces 
to 3.72 A, and from a Ti 3+ to a Ti 4+ cation increases to 
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4.00 A, with respect to bulk structure. 

This type of defected STO becomes insulating upon 
removal of three electrons. The remaing three electrons 
are accomodated in the same number of in-gap localized 
states, as shown in the band diagram of Figure |4j;l, and 
the corresponding PDoS of Figure [5] The localized state 
with higher energy is composed of 4s, 4p z , and 3d z 2- 
orbitals of the Ti 2+ cation, whereas the state with op¬ 
posite spin sign at ~ -0.5 eV below the latter is a com¬ 
bination of the 4s and 3d z 2-orbitals of the same type of 
cation. A state deeper in energy and composed of the 
same orbitals and spin sign than the later is found in 
resonance with the VB states. The resulting net charge 
distribution in the real space displayed in Figure|4^ shows 
a charge density with opposite spin to the localized mo¬ 
ments that it separates which, in turns, vanishes when 
a forth electron is removed from the supercell. Indeed, 
the former state becomes empty with an energy close 
to the CB states (see Figure §) leaving the defected 
area with net charge density joining both Ti 2+ cations 
with no nodes, as plotted in Figure |4ji. Note that one of 
the consequences of depleting electrons from the defected 
structure is an increasing on the total magnetic moment, 
which evolves from zero up to 2 /tb in the four-fold hole- 
doped configuration. 

CONCLUSION 

In spite of the broad interest in STO as a model per- 
ovskite material, a detailed description of the magnetic 
properties of oxygen-defected STO in bulk STO was lack¬ 
ing. We have provided a systematic study of the mag¬ 
netism arising under different doping conditions of mul¬ 
tiple isolated Ov clustered in various geometries to show 
that the local chemistry of the defect is strongly depen¬ 
dent on the defected material doping level. Localized 
magnetic moments may exhibit a large tunability upon 
modification of the charge carriers concentration. Un¬ 
like mono-Ov and depending on the mOv geometry, per¬ 
manent magnetic moments can be induced even in the 
absence of intrinsic doping. These results are consistent 
with recent experiments on complex oxides heterostruc¬ 
tures where the magnetism in STO is strongly dependent 
of the charge carrier concentration with the observation 
of ferromagnetism in electron depleted regions [Bl. The 
understanding and control of the various doping levels 
of Ov and their relative proximity is clearly critical to 
evaluating the possible role of these centers in the mag¬ 
netic behavior of STO under different doping conditions. 
Correlation effects between defected sites can lead to new 


functionalities in designed TMO-based structures. The 
implications of our results may be generalized to other 
defective perovskite-based systems that may now be ex¬ 
pected to exhibit a similar behavior. 

ACKNOWLEDGMENTS 

We acknowledge the computing resources provided on 
Blues high-performance computing cluster operated by 
the Laboratory Computing Resource Center at Argonne 
National Laboratory. Work at Argonne is supported by 
DOE-BES under Contract No. DE-AC02-06CH11357. 
PG was sponsored by the Laboratory Directed Research 
and Development Program of Oak Ridge National Labo¬ 
ratory, managed by UT-Battelle,LLC,for the US Depart¬ 
ment of Energy. J. Low’s technical support is gratefully 
acknowledged. 


* alejandrolb@gmail.com 

[1] A. Ohtomo and H. Hwang, Nature 427, 423 (2004) 

[2] W. D. Rice, P. Ambwani, M. Bombeck, J. D. Thompson, 
G. Haugstad, C. Leighton, and S. A. Crooker, Nature 
materials 13, 298 (2014) 

[3] Ultramicroscopy 86, 289 (2001), international Sympo¬ 
sium on Spectroscopy of Materials 

[4] L. Suescun, B. Dabrowski, J. Mais, S. Remsen, J. W. 
Richardson, E. R. Maxey, and J. D. Jorgensen, Chemistry 
of Materials 20, 1636 (2008) 

[5] K. Kamiya, M. Young Yang, S.-G. Park, B. Magyari-Kpe, 
Y. Nishi, M. Niwa, and K. Shiraishi, Applied Physics 
Letters 100, 073502 (2012) 

[6] F. Bi, M. Huang, S. Ryu, H. Lee, C.-W. Bark, C.-B. 
Eom, P. Irvin, and J. Levy, Nature Communications 5, 
5019 (2014) 

[7] A. Lopez-Bezanilla, P. Ganesh, and P. B. Littlewood, 
ArXiv e-prints(2014), arXiv: 1408.3103 [cond-mat.mtrl- 
sci] 

[8] P. E. Blochl, Phys. Rev. B 50, 17953 (1994) 

[9] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77, 3865 (1996) 

[10] G. Kresse and J. Hafner, Phys. Rev. B 48, 13115 (1993) 

[11] G. Kresse and J. Furthmiiller, Phys. Rev. B 54, 11169 
(1996) 

[12] G. Kresse and D. Joubert, Phys. Rev. B 59, 1758 (1999) 

[13] S. L. Dudarev, G. A. Botton, S. Y. Savrasov, C. J. 
Humphreys, and A. P. Sutton, Phys. Rev. B 57, 1505 
(1998) 

[14] M. Kim, G. Duscher, N. D. Browning, K. Sohlberg, S. T. 
Pantelides, and S. J. Pennycook, Phys. Rev. Lett. 86, 
4056 (2001) 

[15] P. R. C. K. V. R. C. P. Ganesh, H. Zhuang and H. Xu 

[16] Z. Hou and K. Terakura, Journal of the Physical Society 
of Japan 79, 114704 (2010) 



